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']P' We study the cosmological evolution of string/M moduli fields T. We use T-duality to fix the potential and 
Mh. show that the superpotential W^ is a function of the duality invariant function j{T) only. If W is given as a finite 
r^ , polynomial of j then the moduli fields do not give an accelerating universe, i.e. they cannothe used as quintessence. 

Furthermore, at T ^ 1 the potential is given by a double exponential potential V ~ e""*^ leading to a fast 

decaying behaviour at large times. For moduli potentials with a finite v.e.v. of T the energy density redshift is 
;h ' model dependent but if T has a finite mass, ttt- < oo, then the moduli energy density redshifts faster or equal to 

matter. Only if the moduli mass is infinite can the moduli energy density dominate the universe independently of 

the initial conditions. 
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Can string/M moduli fields parametrize a slow rolling cosmological constant? Models with a cosmological 
constant term, intended as a constant vacuum contribution or as a slowly decaying scalar field, have 
recently received considerable attention for several reasons, both theoretical and observational. 

On the observational side, we find direct evidence in recent works on spectral and photometric observa- 
tions on type la supernova |l| that a favour positively accelerating universe. Statistical fits to several 
independent astrophysical constraints support these results [^. See however Q for a different explana- 
tion to these observations. More indirect evidence comes from the observational support for a low matter 
density universe from X-ray mass estimations in clusters 0], H- 

In these works, if the nucleosynthesis limits on the baryonic mass are to be respected, the total matter 
that clusters gravitationally is limited to ^0.3. In such a case, a cosmological term would reconcile the 
low dynamical estimates of the mean mass density with total critical density suggested by inflation and 
the flatness problem. 

The requirements of structure formation models also suggest a cosmological constant term. Simulations 
of structure formation profit from the presence of matter that resists gravitational collapse and A CDM 
models provide a better fit to the observed power spectrum of galaxy clustering than does the standard 
CDM model g, Q. 

From the theoretical point of view, we have to face the possible conflict between the age of the universe 
in the standard Einstein-de Sitter model and the age of the oldest stars, in globular clusters. 

Estimates of the Hubble expansion parameter from a variety of methods seem to point to Hq ~ 70 it 10 
km/s/Mpc (a recent review can be found in ||8[; see e.g. |^ for specific projects), leading to an expansion 
age of i(7 ~ 9 lb 1 Gyr for a spatially flat universe with null cosmological constant. On another hand, 
the age of globular clusters have been estimated in the range ~ 13 — 15 Gyr |jl^; although revised 
determinations based on the Hypparcos distance scale are lower by approximately 2 Gyr |11]. 



Many models with a scalar field playing the role of a decaying cosmological constant "quintessence" have 
been proposed up to now. Some of them are specific models motivated by physical considerations but 
most of them are phenomenological proposals for the desired energy density redshift |12|-[p!^]. 



As a first step in the study of these models, the age of the Universe is calculated for several redshift 
laws of the energy density that resides in the dynamical scalar field |13|, |17], [p!^ ]. Observational con- 
sequences of an evolving A component decaying to matter and/or radiation have been studied in |19|, 



[20 1, obtaining severe constraints on such models. Quintessence [^if effects on the cosmic microwave 



background anisotropy are analysed in |21] and ]2^ and phenomenological difficulties of quintessence 



have been studied in [23|. Constraints on the equation of state of a quintessence-like component have 



been placed from observational data [24|. Recently, a potential for a cosmologically successful decaying 



A term has been constructed in [25|. 



Regarding string/M theory, we have not been able to determine in a unique way the vacuum structure 
of string/M theory |^^, since there are many consistent and physically inequivalent solutions. However, 
recent progress has been made in understanding and connecting many of these vacua using duality 
symmetry and it is hoped that a better understanding of string/M theory will yield to a unique low 
energy limit. 



Even though there are a great number of inequivalent vacua, all these solutions share some general 
properties and the existence of moduli fields is one of them. Our interest in studying these fields is 
twofold since the phenomenological consequences of these fields could be relevant for cosmology and on 
the other hand they could restrict the possible string vacua and in an extreme case rule out string theory 
as a physical theory. 

The cosmological evolution of scalar fields with exponential potentials has been extensively studied in 
the literature pO] , p^ , p^ . These scalar potentials have been motivated by the compactification of the 



extra dimensions and have been used to describe string/M moduli |3C]-[p3]. Using T-duality invariance, 
a valid symmetry in string theory [^I| and M theory |^^, we will show that the string moduli potential 
cannot be approximated by an exponential potential if the superpotential has a finite expansion on the 
modular invariant function j{T) (see below). The resulting cosmological evolution of string moduli will 
then greatly differ from an exponential potential one. 

To begin, we will present the Einstein equations for a flat universe with a Friedmann-Robertson- Walker 
metric in the presence of a barotropic fluid (matter or radiation) and a scalar filed (i.e. the moduli 
field) with arbitrary scalar potential. The cosmological evolution of the moduli field T is determined by 
the asymptotic behaviour of the quantity A = —Vr/V, Vr = dV/dT. We will then proceed to rederive 
and study the non-perturbative scalar potential for the moduli field T and we will obtain the attractor 
solutions. Finally, we discuss the influence of this decaying scalar field on the universe behaviour and we 
present our conclusions. 

Our starting point is a universe field with a barotropic energy density, which can be either matter or 
radiation, and the energy density of the moduli field. The barotropic fluid is described by an energy 
density p^ and a pression p^ with a standard equation of state p^ = (7^ — 1)^^, where 7-^ = 1 for matter 
and 7^ = 4/3 for radiation. We do not make any hypothesis on which energy density dominates, that 
of the barotrpoic fluid or that of the moduli field. For a canonically normalized scalar field (j) with a 
selfinteraction given in terms of the scalar potential V{(f)) but with gravitational interaction with all 
other fields (as moduli fields) the equations to be solved, in a spatially flat Friedmann-Robertson- Walker 
(FRW) Universe, are given by 

p = -3H{p + p) (1) 

where H is the Hubble parameter, T^(</>) is the scalar field potential and we have taken SttG = 1. It is 
useful to make a change of variables [^ x = (j)/\/QH, y = \/V /^/^H and eqs.(|l|) become 



XN 



VN 



-3x + ^j^\y^ + ^x[2x^ + 77(1 - x^ - y^)] 

-\llxxy + ^y[2x^+j^il-x^-y^)] (2) 

Hn = -^H[j^{l-x^-y^) + 2x^] 

where N is the logarithm of the scale factor a, N = ln{a), Jn = df/dN for f = x,y, H and A(A^) = 
—V(j,/V. Notice that all model dependence in eqs.(y) is through the quantities A(iV) and the constant 



parameter 7^. Eqs.(§) must be supplemented by the Friedmann or constraint equation for a flat universe 
^ri + x^ + y^ = 1 and they are valid for any scalar potential as long as the interaction between the 
scalar field and matter or radiation is gravitational only. This means that it is possible to separate the 
energy and pression densities into contributions from each component, i.e. p = p-y + P(f, and p = p-y +P0, 
where p^ (p^) is the energy density (pression) of the scalar field. We do not assume any equation of 
state for the scalar field. This is indeed necessary since one cannot fix the equation of state and the 
potential independently. For arbitrary potentials the equation of state for the scalar field p^ = {'^^ — l)P(j, 
is determined once p^,p^ have been obtained. Alternatively we can solve for x,y using eqs.(|2|) and the 
quantity 'y^ = [p^ +P(i>)/ P<i> = 2x^/(x^ + y^) is, in general, time or scale dependent. 



A complete analysis of the solutions of eqs.(y) has been obtained |15] and the values of x, y at late times 
depend on the asymptotic behaviour of A. If A ^ then x — > 0, y — > 1 and H -^ cte becoming the scalar 
field a "true" cosmological constant. If A — > 00 then depending on whether A oscillates or not we will have 
different asymptotic behaviour. If the v.e.v. is at ^ cxd then A will not oscillate and xocyocA"^ ^0 
and Qfj, — > 0. If, on the other hand, the v.e.v. of (p is finite, A will oscillate and approach infinity. In this 
class of potentials we can expand around the minimum with a leading term V{(j)) = Vo{(f) — 0o)" giving 
an equation of state with 7,^ = 2n/(2 + n) |jT^-[|l^]. If 7</, < 77 ( 70 > 77) then O^ ^ 1 {Q^ — > 0) and 
for 7^ = 77, r2(^ goes to a finite (7^ 0, 1) constant value. We have therefore seen that asymptotic limit 
of the quantity A determines the cosmological evolution of the scalar field and the redshift of the energy 
density. 

Let us now rederive the non-perturbative potential for string moduli. These fields arise in string theory 
when the extra dimensions are compactified and they parametrize the geometry or complex structure 
of the compactified manifold. For string theory compactified on an orbifold there are at least 3 moduli 
T that represent the radii of the manifold {T = R^ + iB with R the radius of compactification and 
B an axion field). Moduli fields do not have any gauge interaction, they interact only gravitationally 
with all other fields and as long as supersymmetry is preserved they have a completely fiat potentials, i.e. 
V{T) = 0. However, once supersymmetry is broken they acquire a steep potential that can be determined 
using T-duality invariance. T-duality has been proven to be present in low energy string theories [31| 



aaaand in M theory |34|. This symmetry ensures that the contribution of the Kaluza-Klein modes, which 
appear because the fundamental object is a string and not a point particle, are taken into account in the 
4-D theory. The T-duality symmetry is given by an SL{2, Z) group where the modulus field T transforms 
as r — > T' = {aT — ib)/{icT + d) with ad — be = I and a,b,c,d e Z |3^. This symmetry gives the famous 
T — > 1/T (i.e. R — > l/R) invariance and it has two duality invariant points T = 1 and T = 6 = e*'^'^ 
plus SL{2, Z) related points. Due to duality invariance we will only need to consider the parameter space 
with T > 1. 

Another important symmetry of the 4-D string theory is supersymmetry. Global supersymmetry is given 
in terms of two functions, the Kahler potential K(T,T), an analytic real function, and the superpotential 
W{T) an holomorphic function that depends only on T (not on T). The Kahler potential gives the kinetic 
term 

Ck = K:f{d^T){d^f) (3) 

with K^ = QrpjSp and the scalar potential is 

V = (K-i)^lWrP (4) 



with Wt = ^^- In order to have a duahty invariant model K and W must be duahty invariant in global 
susy. For K, W duality invariant, the derivatives of K and W are not duality invariant but the scalar 
potential as given in eq.(^) is indeed invariant. For local susy it is only the combination G = K + log\W\'^ 
that must be duality invariant and the scalar potential is given in terms of the covariant derivative 
Gt = ^ = Kt + -y^ by y = e^ {Gt{K~^)'^G'^ — 3). Notice that duality invariance is more restrictive 
for global susy than for local susy since in the former case K, W must be independently duality invariant. 

In string theory, the kinetic term has been calculated in the large T limit and it is given by a non linear 
sigma model 

Ck = l^{d,T){d^f) (5) 

where T^ = T -\- T and from eq.(^) we have Kj^ = 1/T^. The Kahler potential must then he K = 
—log{T + T) + log{f{T)) + log{f{T)). In order to have K duality invariant the function /(T) must 
transform as modular function with weight m = 1 (under a duality transformation a modular function 
transforms as / — > {ic + d)™"/ with m the modular weight) to cancel the transformation of Tr- If we 
do not wish to introduce any spurious zeros or poles in the fundamental domain then / = r](T)'^, where 
r]{T) is the Dedekind eta function defined by r](T) = g^'^^ TlnCl ~ 9")' 1 = e~^ with modular weight 
m = 1/2. 

The superpotential must be a duality invariant function and it is therefore given in terms of the absolute 
invariant function j (T) , since all invariant functions can be expressed as a rational function of j (T) (see 
appendix for definitions of modular functions and useful relations between them but for completeness see 



[35 1 ). All string vacua will then have a Kahler potential and superpotential of the form 

K{T,f) = -log{Tr\v{T)\'), W{T) = W{j{T)) (6) 

respectively. Notice that W depends on T only through j{T). This choice of K and W is valid for global 
and local susy although it is more common to write in local susy Kl = —log{Tr) and a superpotential 
Wl = W ri~'^ with G = Kl + loglWi]"^ = K + log\W\'^. In this later case neither Kl nor Wl are duality 
invariant but G is indeed invariant. 

At large T one has j{T) ~ e^'^'^ and for W = j{T)~"' one obtains a superpotential W ~ g-SawT ^^^ ^ 
scalar potential 

V ~ T^ e-2^''^'- (7) 

leading to the well known exponential potential for moduli fields. Exponential potentials have been 
extensively studied in cosmology and lead to a nonvanishing contribution of the energy density at late 



times and to other interesting properties |Q,|Q,[|2^. However, we would like to point out that for a 



finite expansion of W on j, string moduli fields do not have exponential potentials because their kinetic 
term is non-canonically normalized (c.f. eq.(^)). In fact, taking T real, the canonically normalized field 
(j) = log{T) / \/2 the potential in eq.(|^ is V{(f)) ~ ^^'^'i> g-27rae ^ -^^ j|. jg ^ fiouble exponential potential. 

Different string vacua will have different superpotentials 1^(^(7")), however, we can study them in general 
using the properties of the invariant function j(T) and the physical condition that the potential V should 
vanish at its minimum. This condition is necessary if we want the moduli potential to represent a slow 
varying cosmological constant. If the potential does not vanish at its minimum then we would need to 
explain the smallness of the vacuum energy and we would introduce a fine tuning problem. We therefore 



demand the scalar potential and its first derivative to vanish at the minimum, i.e. V\min = V'\min = 0. 
These conditions require W = as seen from eq.(Q), where from now on the prime denotes derivative 
with respect to T. What are the zeros of W'? We will distinguish two different kinds of zeros: at finite 
T = Tq and at T = oo. The function W' will always vanish at the dual invariant points T = 1, 6 = e^'^'^, 
as any modular invariant function, but it may also vanish at a model dependent point Tq. 

Let us now proceed to study the cosmological evolution of the moduli field. We need to solve eqs.@. Since 
all model dependence in eqs.(0) is given through A = —V'/V its asymptotic behaviour will determine the 



cosmological evolution of the moduli field at late times |15|. We will show that in all cases A will tend to 
infinity for T canonically normalized if the superpotential W is given as a finite expansion on j. Here we 
will mainly consider this kind of superpotentials. For finite values of T, the field will oscillate around its 
v.e.v. and A will also oscillate while approaching infinity. In this class of models the redshift of energy 



density of the moduli field is model dependent ||T4],[15|. On the other hand, for T — > cxo then A — > oo 



without oscillating and the energy density of the moduli field will redshift faster than the barotropic fiuid 
for all models. For the potential given in (Q), we have 

where we have used that Vl^ is a function of j(T) only and Wj = dW/dj. The absolute invariant function 
j{T) has its only zero (triple zero) at T = 6 plus SL{2,Z) related points and it has a pole at the cusp 
point T = oo and by duality at T = 0. The derivative of j, j{T)' , vanishes only at the dual invariant 
points T = 1, (5 and it has a pole at T = oo and T = 0. 

In order to determine the asymptotic value of A at late times we will study the limiting behaviour oi j" / j' 
and j'Wjj/Wj. From the functional dependence of j on T we obtain the following limits (see appendix). 
For T ^ 1 we have j" /f -^ — oo, for T — > 5 we find j"/i' ~^ oo and for a model dependent point 
T = To(/ 1,5) then the value of j" /j' is model dependent but will be, in general, different than zero or 
infinity. Finally, for T ^ oo we have a constant value, j" /f — > 27r. 

Now, we consider the last quantity in eq.(P). Without loss of generality if < T > is finite, we can expand 
W around the zeros of W' that we wish to analyse, i.e. W = T,nan{j{T) — jo)" with j{Tq) = jo and 
V{To) = W'{To) = 0. Since in the asymptotic limit T ^ Tq and j{T) -^ jo the leading term of Wj and 
Wjj is the one with the smallest exponent, say no, giving Wjj/Wj ~ ao{no — l)/(j — jo)- For Tq ^ 1,6 
but finite we have j' ^ 0, and j'Wjj/Wj ^ oo in the limit T -^ Tq. If the minimum is at T = 5, then 
j(ro) = jo = and since j has a triple zero while j' only a double zero j'Wjj/Wj — > ao(no — l)j" /j' 
also approaches infinity while for Tq = 1 we have two possibilities. Since W'(T = 1) = j'{T = 1) = 
for any value of jo we could take jo = and then j'Wjj/Wj -^ 0. On the other hand, we could take 
Jo = i(l) = 1728 in which case j'Wjj/Wj -^ oo. Finally, for T — > oo we find j'Wjj/Wj = ao{no — 1)2tt 
which is the only case where it remains finite. 

We have found, so far, that A ^ oo for finite values of T (i.e. T = 1,5, To) and it tends to a constant 
A -^ — 27r(ao(no — 1) — 1) for T -^ oo. However, the moduli field T has non canonically normalized kinetic 
terms and eqs.(|l|) or (^) are valid only for canonically normalized fields. The canonically normalized 
field is ^ = /n(T)/v2 (we are taking T to be a real field since its imaginary part is an axion field). In 
terms of (/) the relevant cosmological parameter is Xc.n. = —{dV/d(j))/V = —\/2T\ giving Xc.n. ^ oo in 
all cases, i.e. for T -^ 1,S,Tq,oo. Besides the asymptotic limit it is also the oscillating behaviour of A 



that is important in determining the cosmological evolution of the moduh field. In table 1 we show the 
behaviour of several cosmological relevant quantities: Xc.n.,7T,^T = Pt/p and a for the different v.e.v. 
of T. The acceleration parameter a is defined as |15] 

^ p + 3p ^ 37-2 

" - (37^ - 2)p 37^-2 ^ ^ 

with J = {p + p)/p- If a = 1 then the acceleration of the universe is the same as that of the barotropic 
fluid and any deviation of a from one implies a different cosmological behaviour of the universe due to the 
contribution of the scalar field. A positive accelerating universe requires a negative a while for < a < 1 
the acceleration of the universe is negative (deceleration) but smaller than that of the barotropic fluid. 
For a > 1 the deceleration is larger than for the barotropic fluid. In terms of the standard deceleration 
parameter q = — 1§ one has a = 37-^ or in terms of x,y one finds a = 1 — 377^(2/^ 7^2;^) = 



77 -7r 



1-3^^377-2- 

Our analysis so far has considered finite expansion of the superpotential only. However, a quite interesting 
case arises for a superpotentialQ W = Log[j{T)]'^. In this case, the leading term in j and j' for T ;::^ 1 in 
the scalar potential V = a^T^lf /j\'^\Log{j)\'^^"-^^^ cancels out and V ~ T"^"- = e"^^'^. This superpotential 
is the only one which allows for an exponential potential for the moduli field and it has Xc.n. = — a2v2. 
The vanishing of the potential at infinity requires a < 0. 

When the minimum of the potential is at T = 00, A approaches infinity and does not oscillate. In this 
case the energy density of the moduli tends to zero and x, y ^ in a few e-folds. This result excludes 
string moduli field as a slow varying cosmological constant when the moduli field tends to infinity. It 
is common in the literature to use exponential potentials to parameterize moduli fields in the T ^ 1 
region, however, we have shown that for finite expansion of W on j, the correct asymptotic behaviour 
of moduli fields is to have a double exponential potential {V = Vqc""'^ ) giving a completely different 
cosmological evolution for the moduli field. In fig.(l) we show the evolution of different cosmological 
parameters for a superpotential W = j{T)^^ which gives a scalar potential V = T^|//j~^p. In the 
region T ^> 1 one has V ~ T^e"'^'^ , T = e^^*^. We can see that the behaviour of x,y is to rapidly 
approach zero with A constant. This is the scaling regime. Afterwards, A increases and tends to infinity 

3" 77 ,, ^ . / 3(2-7^)7. " 



while x,y after a subsequent increase end up going to zero with 2; — > a/| -^, y — > V ~2^^^ Ell' "^^^ 
behaviour of the moduli field with the complete duality invariant potnetial is very similar to the double 
exponential case shown in [^] . The evolution of the Hubble parameter H approaches Hm (i.e. for matter 
fields) after a few e-folds while the accelerator parameter a goes to one. The acceleration of the universe 
is then not affected by the moduli field. The same conclusion holds for the expansion rate of the universe 
given by 7 = p/{p + p) = {a - 2)/3. 

If the v.e.v. of T is finite the cosmological evolution of the moduli field is model dependent. We have 
seen that in all these cases A — > 00 with an oscillating behaviour. This is of course due to the oscillating 
behaviour of the field around the minimum of the scalar potential. Without loss of generality we can 
expand the scalar potential around the v.e.v. of the canonically normalized scalar field (p and keeping 
only the leading term one has ^ ~ Vo((/> — 0o)"- This kind of potentials have an average equation of 
state parameter 7^ — > 2n/(2 + n) [^,jl^]. The condition to have a finite mass (i.e. ^^ = m? < 00) 



^A complete analysis will be presented elsewhere 
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Figure 1: We show different cosmological parameters for a potential V — cT^\j' /j~'^\'^ with T ^ oo at late times. The 
initial conditions are xo = 0.1, j/o = 0.5 and Hq = 1. In fig. 2a we show the evolution of Qt- In fig. 2b we see the behaviour 
of X (dot line) and y (solid line) . In fig.2c we plot the Hubble parameter (dash-dot line) compared to a matter dominiated 
universe Hm (dot lino). In fig. 2d we have Ac.n.. Finally we show in fig.2e and 2f the 7c. n. a parameters respectively. 



requires n > 2 which imphes that 70 > 1. This means that the redshift of the moduh field will be faster 
or at best equal to that of matter fields and therefore it will not parameterize a slow rolling cosmological 



constant. However, if we don't impose a finite moduli mass then the condition V\ 



V'\ 



only 



imposes the constraint n > 1. For 1 < n < 2 we have 2/3 < 7,^ < 1 giving a moduli redshift slower than 
for matter fields. So, even though a slow varying cosmological constant can be obtained and it could help 
to solve the universe age problem, the condition n > 1 implies that it will not positively accelerate the 
universe (i.e. a < 0) as required by the recent supernova results. For n < the moduli field will tend 
to infinity and the correct scalar potential is given by the double exponential potential (see above). In 
terms of a superpotential W = {j{T) — jo)"" and the scalar potential V ~ Vo(' 
between no and n is: n = 2 for for T = 1 independently of the value of no, n 
and n = 2(no — 1) for a model dependent value To. The difference in the relationship between n and 
no is due to the different multiplicty of the zero around the value of To that we are expanding. If the 
superpotential has a finite expansion on j then the moduli fields will not have a scalar potential with 



(/>o)" the relationship 
2(3no - 1) for Tq = S 



negative powers of T, i.e. V ~ T , 6 < 0, since the leading term in j and / is given by 1/q 
at large T any dependence in V on j or j' will give the main contribution. 
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Table 1. We show the asymptotic behaviour of Ac.n., 7t, a and JIt for different v.e.v. of T and a superpotential with leading 
term W = {j{T) —jo)^° and we have taken 7-y = 1 in a. In the last line we have taken a superpotential W = Log[j{t)]°', a < 
and the first quantity in 7t, a, fir corresponds to a^ > 37^/8 while the second one for a^ < 3/4 respectively. 

We give two examples of oscillating moduli fields. In fig. 2 we show the evolution of different cosmological 
parameters for a superpotential W = j{T)^ and V = cT^\2j'j\'^, with^c = 10~^^ and initial conditions 
xq = 0.1, yo = 0.5, i/o = 1- The scalar potential vanishes at the dual invariant points and we study the 
behaviour around To = 1. Since in this case V cc {(p — 1)^, i.e. n = 2, the moduli energy density should 
redshift as matter in the asymptotic region and as can be seen from fig. 2a. We also see in fig. 2b that 
X, y approach zero quite rapidly and they start oscillating when the scaling period is over, as can be seen 
by comparison with fig.2e. The acceleration parameter tends to one after a few e-folds and 7t ends up 
oscillating with an average < 7^ >= 1. 

As our second example for an oscillating field we show, in fig. 3, the cosmological evolution for a super- 
potential W = {j - jo)^/^^ with a scalar potential V = c{8/15fT^\j'\^\{j - io)^'^/^^ with jo = j{T = 
2), c = 10~^ and initial conditions xq = 0.1, yo = 0.5, Hq = 1. The scalar potential vanishes at the 
dual invariant points and at Tq = 2 and we study the behaviour around Tq = 2. Since in this case 
V oa {4> — 1)^'^, i.e. n = 6/5, the equation of state parameter is < 7t >= 3/4 and the moduli energy 
density must redshift slower than matter. This can be seen in fig. 3a. Notice in fig.3e and fig.3f that the 
average jt and a is different than one leading to a different behaviour of the universe with the presence 

We took a small c to allow the initial value of Ti not to be too close to To = 1. For c = 1 one has Ti = 1.007 while for 
c = 10"^^ one has Ti = 2.04. 
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Figure 2: We show the evolution of different cosmological parameters for a potential V — cT^\2j'j\^ with T — + 1 at late 
times. The initial conditions are xo = 0.1, t/o = 0.5 and Ho = 1. In fig. la we show the evolution of fir- In fig. lb we see the 
behaviour of x (dot line) and y (solid line). In fig.lc we plot the Hubble parameter (dash-dot line) compared to Hm (dot 
line) and the a parameter (solid line). In fig. Id we have jt- Finally we show in fig.le and If the Xc.n. and the evolution of T. 



of the moduli field. This can also be observed by the discrepancy between H and H^ in fig. 3c. However, 
it is important to remember that in this case 1^'| = but V"\ — > oo giving an infinite mass for the moduli 
at < T >= 2. 

To conclude, we have studied the cosmological evolution of string moduli fields in a completely general 
way. We have only assumed the well established results of T-duality and the sigma model kinetic term for 
the moduli field (at large T). The scalar potential is model dependent but duality invariance constraints 
the potential severely. If the superpotential is given as a finite expansion on j then, in the region of 
large T, the scalar potential has a double exponential form leading to an energy density redshift of the 
modular field faster than the barotropic fluid. This result excludes this kind of potentials as a slow 
varying cosmological constant. Furthermore, for a finite v.e.v. of the moduli, if its mass is finite (i.e. 
7n(To) < oo), then the redshift of the moduli will be no slower than matter fields and will not be able 
to contribute to solve the age crises nor the positively accelerating universe. However, if we allow the 
moduli mass to be infinite then the redshift will be slower than for matter fields and it will give the 
main contribution to the energy density of the universe but it will not lead to a positively accelerating 
universe. We have therefore shown that moduli field can not parameterize a slow rolling cosmological 
constant with a positively accelerating universe. 

I am greatful to G. Piccinelli for many useful discussions. This research was supported in part by 
CONACYT project 32415-E and by DGAPA, UNAM, project IN-103997. 



Appendix 

In this appendix we give some useful relations between modular functions |B5|. A modular function / with weight m 
transforms under duality SL{2,Z) as /(T) ^ {icT + d)"" f {T) for T ^ {aT ~ ib) / {icT + d) , ad - be = 1, a,b,c,deZ. This 
symmetry has two dual invariant points T = 1,5 = g*'^/^ plus SL{2,Z) related points. A complete set of independent 
modular (holomorphic) functions is given by the Eisenstein series, with modular weight 2k, 

-2k 



G2k(T) = E„i,„2jz(mir + 712) 

.(T) = 2C(2fc)(l-i^E^=ia2fc-i(n)q") (10) 



G2k( 

where q = e~'^^'^ , C, is the Riemann's zeta function, B2k are the Bernoulli numbers and ap{n) is the p-powers of all divisors 
of n. The Eisenstein functions (for fc > 1) do not have any poles and their zeros are at T = 1 for k odd and at T = 5 
for k 7^ OmodS. For k — 1 the order of summation must be specified and one gets either G2{T), which is holomorphic 
but transforms as a connection G2 — > {icT + d)^G2 — 2Tvic{icT + d), or G2{T,T) — G2(T) — iv/ReT, which transforms as 
a modular function with weight two but it is not holomorphic. G2 does not have any poles and its zeros are at the dual 
invariant points T — 1, S — e"^' . 

The derivatives of G2k for k > 1 and k — 1 are given by 

k 
G'2k{T) = a{k)G2k+2 G2G2k 

TV 

G',(T) = ^G.-i-G^ (11) 

with a(k) = ■KkC,{2k) /ZC,{2k + 2). The derivative of any modular invariant function vanishes at the self dual points. We can 
express the Dedekind Eta function ri(T), the absolute invariant function j{T) and their derivatives in terms of the Eisenstein 
functions 

7,(T) = Ai/2^ = g^/2*n„(l-g"), ^'(T) = -^r,G2 (12) 
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Figure 3: We show the evolution of different cosmological parameters for a potential V = cTr\j'\^\{j — io)|~^ with 
Jo = j(2) and T ^ 2 at late times. The initial conditions are xo ~ 0.1, yo = 0.5 and Ho = 1. In fig. 3a we show the evolution 
of JIt. In fig. 3b we see the behaviour of x (dot line) and y (solid line). In fig. 3c we plot the Hubble parameter (solid line) 
compared to a matter dominated case Hm (dot line) while in fig. 2d we have Xc.n.. Finally we show in fig.le and If the 7t 
parameter and the acceleration parameter a respectively. 
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with the discriminant 



JiT) = ■ 
A(T) 



A' 



675 



j'iT) 



21 .Ge 
■K^ Ga 



{2QGI - 49G6) 



(13) 



(14) 



256 7ri2 

The 77 function has a zero at T = cxa and its dual invariant point T = Q and it has no poles. The invariant function j{T) has 
a triple zero at T = 5 and it has a pole &t T = 00. Using eqs.(|ll|) and (113) we can get some useful relationships 



G6(^) = ^(,-1728)^/^ 



7 2 
G4(T) = — Ge G2G4 

TT TT 

10 3 

G(j(T) = — Gs G2G6 

TT TT 



(15) 



For finite T, G4 has its only zero at T = 5 while Gs at T = 1 and they do not have any poles. Using the properties of the 
invariant function j and of the Eisenstein function G2k one finds 



The expansion around T ^ 1 gives 



f ^ 21 Gs ^ 10 Gs _ 7 Gs _ J_^^ 

j' TT Gi TV Ge TT G4 TT 



G2{T) = y(l-24Q-72g2 



^(l + 240g + 2160g'+, 

— -(l-504q-16632g^ 
945 



G4(T) 

GeiT) 

J{T) = ■ 
y 

ri{T) = q^^'"'{l + q + 2q^+3q^ + 
To summerize, the values of G2, G4, Ge,ri, j, j' ,j" at T — 1,5,00 are 



744 + 196884g + . 



T 


G2(r) 


G4(r) 


G6(r) 


V{T) 


JiT) 


j'(r) 


J"{T) 


1 


TT 


3.15 TT-* 





0.768^2^ 


1728 





49655 


S = e*^/« 


27r/V3 





5.86 TT*^ 


(0.79 -i 0.10)^2^ 


(triple) 


(double) 





CXJ 


7rV3 


7rV45 


27rV945 


q'/'^ - 


1/g -^ cxj 


27r/g -^ 00 


47r^/g -* cxj 



(16) 



(17) 
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